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Abstract 

I discuss effective field theories for heavy bound systems, particularly bound systems involv- 
ing two heavy quarks. The emphasis is on the relevant concepts and on interesting physical 
applications and results. 

1 NON-RELATIVISTIC BOUND SYSTEMS 

In nature there are many particle bound systems for which the relative velocity v of the particle in 
the system is small. Typical examples in the domain of electromagnetic interactions are positro- 
nium (e + e~) and muonium (e~/u + ), their counterpart in the strong interaction domain are heavy 
quarkonia (ti, bb, be, cc) and somehow in the middle lie systems like hydrogen, hydrogenoid atoms, 
pionium {n + n~). All such systems are non-relativistic bound systems. 

At the beginning of this paper, I will use the example of positronium to show what are the 
typical (technical) problems inherent to a non-relativistic bound state calculation, even in a pure 
perturbative situation and I will relate such problems to the existence of several physical scales. 
Then, I will discuss the further complications that arise in bound state calculations inside a strongly 
coupled theory like QCD. For the rest of the paper, I will introduce non-relativistic Effective Field 
Theories (EFT) and I will explain how EFT greatly simplify the bound state calculations, both 
technically and conceptually, allowing us to obtain systematically interesting and new physical 
results. 



2 AN EXAMPLE of bound state dynamics: POSITRONIUM 

To understand what are the peculiarities of the bound state interaction, let us first consider positro- 
nium. Here, the interaction causing the binding is the electromagnetic interaction: QED fully 
describes this system and the coupling constant is the fine structure constant a, which is small and 
completely under control in the region of physical interest. The positronium energy levels are given 
by the poles of the four-point positron-electron Green function G = (OlV'i V 7 2' i / ; 2' ? / ; i |0) ? which in turn is 
given as a formal series in terms of Feynman diagrams and thus in terms of a. An integral equation 
can be written, called Bethe-Salpeter equation, whose solution is the four-point Green function G 
and whose kernel K is the subset of two-particle irreducible diagrams. Appropriate techniques allow 
us to write down formally the energy levels as a series of contributions involving kernel insertions on 
the zeroth-order Green function averaged on the zeroth order wave functions (for a general review 
see [0], for an explicit calculation see e.g. Q). 

There is an important difference between the calculation of a scattering amplitude and the 
calculation of bound state wave functions and energy levels. Both are given in terms of Green 
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functions projected on initial and final states, but in the first case the initial and final states are 
on shell, i.e. the wave functions describing initial and final particles are free, while in the second 
case the initial and final wave functions are the bound state ones and thus bear a dependence on 
a. This last fact produces that from the expansion of the energy levels it is not trivial to select the 
set of diagrams contributing at any given order in a. In fact, being the wave function a-dependent, 
the number of vertices in a diagram do not allow to trace back the order in a of the contribution of 
the diagram. Precisely, it happens that: the contribution of each graph is a series in a; the leading 
order in a does not follow from the number of vertices in the graphs. For example, diagrams 
differing only for the number of ladder photons that they contain, contribute all at the same leading 
order in a to the energy levels, plus subleading contributions. Therefore, it is necessary to resum 
an infinite series of diagrams to complete an order in a in the energy levels. In other words, the 
bound state calculation are 'non-perturbative' in the binding interaction. Moreover, the contribution 
of a diagram depends strongly on the gauge. In particular, spurious terms can be generated and 
canceled only by subsequent contributions. In fact, while in a scattering calculation gauge invariance 
is manifest term by term at any order of the expansion in a, in a bound state calculation one has 
to look for particular subsets of diagrams and prove that they are gauge invariant. This turns out 
to be quite difficult in practice. 

The same fact can be explained in the following alternative way. Let us consider the calculation of 
the perturbative corrections to the electron anomalous magnetic moment. This is a typical scattering 
amplitude calculation and in the integrals related to the actual evaluation of the Feynman diagrams 
we have only one relevant physical scale: the mass m of the electron. On the other hand, let us 
consider the calculation of the positronium energy levels. This is a typical bound state calculation 
and in the integrals related to the actual evaluation of the Feynman diagrams we have three relevant 
physical scales: the mass m of the electron, the relative momentum p ~ mv and the bound state 
energy E ~ mv 2 . For positronium v ~ a <C 1 and thus such scales are quite different and get 
entangled in the calculations causing the Feynman diagrams to contribute in a nontrivial way to 
the perturbative expansion in a, as explained above 

This is what happens in the QED calculation of the bound state energy levels. Alternatively, 
one can take explicitly advantage of the non-relativistic nature of the system and start directly 
from the non-relativistic reduction of the QED bound state problem. In this case the zeroth-order 
problem is the Schrodinger equation 

f- — +V(r) -E)G'(t,t',E) =5 3 (r-r'), (1) 
V m ' 

where V(r) is the Coulomb potential V = — j and G' is now the non-relativistic Green function. 
The Coulomb potential term corresponds to the resummation of all the ladder photon contributions. 
In other words, Eq. (|l]) provides a resummation of the most singular terms ~ (^) n in the Feynman 
diagrams series. Therefore, this appears to be the suitable starting point to study a non-relativistic 
bound state like positronium. However, in order to proceed further and systematically calculate 
perturbative and relativistic corrections, the non-relativistic reduction needs to be under control. 
A number of problems emerge in any kind of naive non-relativistic reduction like: 

• Considering retardation (or nonpotential) effects. Problems typically arise on trying to go 
beyond the leading order (|l|). Retardation effects, related to low energy photons (or gluons 
in QCD), appear at some point in the reduction. Any precision perturbative calculation 
has therefore to deal with them. However, no systematic procedure was developed before 
the effective field theory approach. Such an issue becomes particularly relevant in QCD 
where nonperturbative contributions appear also as nonpotential effects (due to very low 
energy gluons) ||. The lack of a systematic and clear approach to evaluate potential and 
nonpotential effects, led in the past to many inconsistencies and contradictory statements 
(about the existence or non-existence of the qq potential) in the frame of the naive non- 
relativistic reduction. 



• Considering relativistic corrections. Problems typically arise on considering the contributions 
of the relativistic corrections in the 1/m expansion to the energy levels or to the wave function. 
In this case ultraviolet (UV) divergences may arise at some order of the (quantum mechanical) 
perturbative calculation. The UV divergences reflect the fact that the naive non-relativistic 
limit is just an approximation of the field theory, valid for small moment of order mv or 
smaller. 

• Calculating contributions in quantum mechanical perturbation theory. As soon as subsequent 
contributions in the 1/m expansion are obtained, a question arise about the treatment of 
such terms. Have we to include as many terms as possible inside the Schrbdinger equation or 
have we to treat them perturbatively? Such issues are often encountered in the literature for 
example in connection to the spin-spin term ~ ^z^i ' S2<5 3 (r). Being such a term not bounded 
from below, the delta function is often substituted by a suitable regularization and then used 
inside the Schrodinger equation instead that in (quantum mechanical) perturbation theory. 

In the next Sections I will show how the effective field theory formalism provides a definite solution 
to these problems allowing us to fully exploit the simplifications related to the non-relativistic nature 
of the interaction. 

Here, I would like to point out that the difficulties discussed above belong to bound state calcu- 
lations in a theory fully under control like QED. When we start addressing a strongly interacting 
theory like QCD in the low energy region where the heavy quark bound state lie, new conceptual 
complications arise. They are related to the existence of nonperturbative contributions and to the 
presence of a new physical scale, Aqcd, the scale at which nonperturbative contributions become 
dominant. In this situation, both the Bethe-Salpeter approach || and the naive non-relativistic 
reduction Q are based ad-hoc ansatz and drastic as well as non systematically improvable approx- 
imations and are biased by ambiguities and indeterminations. In the case of heavy quark bound 
states, the effective field theory approach provides us with a systematic and under control approach 
where high energy (perturbative) and low energy (nonperturbative) contributions are disentangled. 
This is what I will discuss in the next Sections. 

More in general we can say that non-relativistic effective field theories give an answer to the 
long-standing problem of how to derive quantum mechanics from field theory. 

3 TAKING ADVANTAGE of the PHYSICAL SCALES 

Let us consider first heavy-light mesons. In these systems there are two relevant physical scales: the 
mass m of the heavy quark, which is very large and in this sense is a 'perturbative' scale, and the scale 
Aqcd which is low and is a 'nonperturbative' scale. Apart from the mass, all the other dynamical 
scales, like the momentum p and the energy E, reduce to Aqcd m heavy-light systems. Heavy 
quark effective theory (HQET) || is the QCD effective field theory that allows us to separate hard 
momenta (p ~ m) and soft momenta {p ~ Aqcd)- Hard effects can be calculated in perturbation 
theory; soft effects are governed by the spin- flavor symmetry that makes the theory predictive. 
Calculations are organized in expansions in a s (m) and Aqcd/ 7 ^- The equivalence between HQET 
and QCD is imposed via the 'matching' procedure. Such procedure has been carried out in many 
(equivalent) ways: by imposing the off shell Green functions to be equivalent to those of QCD 
0], by integrating out the 'antiparticle' degrees of freedom ||; by performing a Foldy-Wouthysen 
transformation in the QCD Lagrangian |J. Precisely, the HQET is organized in a power series in 
the inverse of the mass of the heavy quark. Only two kind of terms contribute: 1) terms containing 
light degrees of freedom (gluons and light quarks) only; 2) terms containing a bilinear in the heavy 
quark field. The size of each term is estimated by assigning the scale Aqcd to whatever is not a 
heavy mass in the Lagrangian. The effective theory thus provides us with an unambiguous power 



counting for each operator and allows us to make precise calculations accurate up to the desired 
order of the expansion. 

The physical situation for bound states formed by TWO heavy quarks is more complicate. We 
have now four scales into the game: m, p ~ mv, E ~ mv 2 and Aqcd- In QCD it is not true 
that v ~ a s always. Such relation holds true only in the strictly perturbative regime defined by the 
condition: mv ^ mv 2 ^> Aqcd- In general v is a nonperturbative quantity, a function of both a s 
and Aqcd: v = v (a s , Aqcd)- However, the only relevant fact is that, being w« 1, the different scales 
turn out to be well separated. The presence of such entangled scales makes perturbative calculations 
very difficult as it was explained before. Moreover, it makes even a numerical calculation beyond 
reach. In fact in a lattice simulation of a heavy quark bound state, one should have a space-time 
grid that is large compared with 1/mv 2 but with a lattice spacing that is small compared with 1/m. 
To simulate e.g. bb states where m/mv 2 ~ 10, one needs lattices as large as 100 4 which are beyond 
our present computing capabilities |Iq| . 

For all these reasons, it is relevant to construct effective field theories which are able to dis- 
entangle these physical scales. One can think about constructing an EFT where v is the small 
nonperturbative expansion parameter and the correction in the expansion are thus labeled by pow- 
ers of v (power counting in v). The first effective field theory of this type that has been introduced 
is non-relativistic QCD (NRQCD) (and NRQED for QED)|o|. NRQCD factorizes the scale m but 
does not deal with the scale hierarchy mv ^> mv 2 . Three scales p, E, Aqcd remain therefore en- 
tangled and even in the perturbative regime, complications arises in bound state calculation due to 
the fact that the power counting becomes not univocal beyond leading order. The true problem is 
that NRQCD does not make explicit the dominance of the potential interaction (in the perturbative 
regime the dominance of the static Coulomb force) and the approximate quantum mechanical nature 
of the system, as has been discussed in Sec. 2. Then, another EFT that enjoys these characteristics 
has been introduced and is called potential Non-Relativistic QCD (pNRQCD). In the next sections 
I will introduce NRQCD and pNRQCD. For the rest of the paper I will mainly present and discuss 
results and predictions of pNRQCD. 



4 NRQCD 

The mass m can be removed from the dynamical scales of the problem using renormalization tech- 
niques. The idea is to introduce an ultraviolet cutoff A of the order of the mass m or less (but much 
larger than any other physical scale). This cutoff explicitly excludes relativistic heavy quarks from 
the theory. On the other hand this is a sensible choice of cutoff since the physics of heavy quark 
mesons is dominated by momenta p ~ mv. Since we are dealing with a quantum field theory, the 
relativistic states we are cutting out have a relevant effect on the low-energy physics. But we can 
compensate for this loss by adding new coefficients and new local interactions to the Lagrangian. 

To leading order in 1/A such interactions are identical to the interactions already present in 
the theory; beyond leading order we have to include nonrenormalizable interactions multiplied by 
some 'matching' coefficients Cj . In principle there are infinite such terms to be included, in practice 
we need only few of them. If we aim at an accuracy of order {p/ K) n we keep in the Lagrangian 
only terms up to and including 0(1/A) n interactions. The coupling m, g, Cj are determined by the 
requirement that the cutoff theory reproduces the results of the full theory up to order (p/A) n . In 
practice is convenient to take A ~ m, to organize the Lagrangian in powers of 1/m making thus 
explicit the non-relativistic nature of the physical systems and to transform the Dirac field so as to 
decouple its upper components from its lower components, i.e. separating the quark field from the 
antiquark field (Foldy-Wouthysen transformation). 

Equivalently, we say that by integrating out, in the sense of the renormalization group, high 
energy degrees of freedom we pass from QCD to NRQCD. The technical procedure to achieve 
this passage is called "matching" and the information about the integrated degrees of freedom is 



contained in the " matching coefficients" . The matching scale separates the degrees of freedom we 
integrate out in the EFT from those that remain dynamical in the EFT. Such scale corresponds 
to the cutoff of the EFT if the regularization procedure is based on cut-off regularization. At each 
matching step the non-analytic behaviour in the scale which is integrated out becomes explicit in 
the matching coefficients. In this case we are integrating out the mass that becomes an explicit 
parameter in the expansion in powers l/m n in the Lagrangian while the dependence in ln(m//n) is 
encoded into the matching coefficients. 

Only recently it has been realized how to perform the matching between QCD and NRQCD 
in dimensional regularization and it has been understood that it is the same matching as in 
HQET (plus four fermion terms)[]. In fact it is true that, since two scales p ~ mv and E ~ mv 2 
(plus Aqcd) exist in NRQCD, while only one scale is left in HQET p ~ Aqcd, E ~ Aqcd, the 
power counting in the two EFTs is completely different. However, we should not be mistaken by 
this fact. The key point here is that in the matching, the power counting of the effective theory 
plays no role at all, the only relevant thing being that all the other dynamical scales (present in the 
effective theory) are much lower than the matching scale. The matching computation in the effective 
theory in dimensional regularization is zero, which greatly simplifies the calculation. Dimensional 
regularization, and the MS scheme is a useful tool that we will always imply in in the following. 

Since the scale of the mass of the heavy quark is perturbative, the scale \i of the matching 
from QCD to NRQCD, mv < /x < m, lies also in the perturbative regime. Then, it is possible 
"to integrate out" the hard scale by comparing on shell amplitudes, expanded order by order in 
1/m and in a s , in QCD and in NRQCD. The difference is encoded into the matching coefficients 
that typically depend non-analytically on the scale which has been integrated out (like ln(^/m)), 
in this case m. We call this type of matching perturbative matching because it is calculated via an 
expansion in a s . 

Up to order 1/m 2 the Lagrangian of NRQCD [1C, 12, 13] reads: 

t / D 2 D 4 S B D E-E D 

^NRQCD = V { hA) + C 2 — + c 4^3 + Cp 9—^- + C D9 ^2 ( 2 ) 

. S-(DxE-ExD)\ , , . , , i 

+tcsg — ^ ip + antiquark terms + terms with light quarks 

4m z J 

_^L/? a J7 a ^^]^p a r) 2 F a ^ u + —nf i, F a F b F c 

m z m z m z m z 

ip and x being respectively the quark and the antiquark field; is the covariant derivative, E and 
B are chromoelectric and chromomagnetic fields, F^ v is the gluon field strength, S is the total spin, 
T a is the color generator, g is the coupling constant, a s = g 2 /4ir. Ci, bi,di are matching coefficients. 
The gluonic part in the third line comes from the (heavy quark) vacuum polarization while the last 
line contains the four quark operators. 

Since the mass m is explicitly removed from the dynamics, a lattice evaluation of a typical heavy 
quark bound system like bottomonium may now take place at lattice spacings larger by a factor 
1/v. This reduces the needed size of the lattice by a factor 1/v 4 , which is approximately a factor 
100 for the IS of bb. 



Actually only the phenomenological use makes the difference between HQET and NRQCD. In particular, at the 
level of the Lagrangian, the power counting is different. For example, terms with two heavy and two light quark 
fields are relevant in heavy-light systems but not in heavy-heavy systems and thus are considered in HQET and not 
in NRQCD. Corrispondently, terms with four heavy quark fields are relevant only for heavy-heavy systems and not 
for heavy-light and thus they are considered in NRQCD but not in HQET. The HQET Lagrangian obtained with 
matching exploiting procedures different from the Foldy-Wouthysen expansion is in fact different from the NRQCD 
Lagrangian but it is related to it via local field redefinitions or via the equations of motions. 
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Figure 1: Matching procedure: on shell green function expanded in a s and 1/m are equated in QCD 
and in NRQCD. The difference is encoded in the matching coefficients. 



The matching coefficients depend on fj, and m and are known in the literature at a different 
level of precision ||i~3| , 14, p!7f| . The [i dependence in the matching coefficients cancels against the 



/i dependence of the operators in the Lagrangian. In particular, 61 fixes the new /^-function ][T^]. 
Let us see in which way. Suppose that in QCD we have nj flavors. Then in NRQCD we have 
y = rif — 1 relativistic flavors. Consequently the MS running coupling constant in NRQCD is 



11 



expected to run according to nj — 1 flavors. In fact in (j2|) the gluonic kinetic term does not have 
the standard normalization but it is multiplied by b\. This can be recovered by a field redefinition 
of the gluon field. Since the remaining gluon field in Lnrqcd are multiplied by g, this is equivalent 
to make in these terms the change 



Q:, 



f-i as rr 1 m A 1 as 



(3) 



At one loop this corresponds to substituting the running coupling constant of nf flavors with the 
running coupling constant of nj — 1 flavors. Therefore, the a s multiplying the gluon fields is 
understood at the scale \i running according to rif — 1 flavors. In all the matching coefficients not 

nf can be substituted with a^ 1 1 (being the difference in higher orders 



multiplying gluon fields, a s 
in a s ) and the scale turns out to be fixed to m 



a s (m). 



The coefficients 62, and d\ . . . are known at 1-loop, |ll], |13J and cf at two loops (in the 

anomalous dimension) WM, 
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(4) 



being 70, 71 the anomalous dimensions, and |15|] 

1 
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q s (m)\ 13C ^ 6 ^' 
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a s (m) \ 2CA/l30 



(5) 



where we have presented the renormalization group (RG) improved results. Notice that the one loop 
correction (finite part) to the chromomagnetic coefficient cf turns out to give a large contribution, 
of order 15% in bb states and of order 30% in cc states. For what concerns the evaluation of 
the spin splittings therefore the matching coefficient is more important than the subsequent 0(v 4 ) 
corrections. 

The effective field theory must have the symmetries of the fundamental theory it is derived from. 
Hence NRQCD/HQET must have the symmetries of QCD, in particular Lorentz invariance. Such 



symmetry is implemented in a nontrivial way as a reparametrization invariance [16]. In practice, thus 
Lorentz invariance is realized through the existence of relations between the matching coefficients, 
precisely: C2 = C4 = 1, cs = 2cf — 1- We stress that : 

• The matching from QCD to NRQCD is perturbative. 

• In NRQCD two dynamical scales are still present, mv and mv 2 . Because of the existence of 
these two scales, the size of each term in (|2|) is not unique. Counting rules have been given fi"9|] 
to estimate the leading size of the matrix elements of the operators in (pi) but still they do 



not have a unique power counting but they also contribute to subleading orders in v. No 
explicit power counting rules to systematically incorporate subleading effects exist, even in 
the perturbative situation. The entangling of the two surviving dynamical scales makes still 
very difficult the evaluation of the Feynman diagrams in any perturbative calculation. 

It is a consequence of the previous observation that the NRQCD Lagrangian contain, mixed, 
both potential and retardation corrections. Indeed both potential (or soft (S)) and ultrasoft 
(US) degrees of freedom remain dynamical in NRQCD, while the scale of the mass (hard) has 
been integrated out. 

The rest mass m has been removed from the quark energies, allowing for much coarser lattice 
than in the Dirac case. The quark and the antiquarks have been decoupled and thus the 
quarks Green function satisfies a Schrodinger like equation 



(iD 



+ 5! + . . . ) G NRQCD (X; xl) = 6 4 {x _ x>) (6) 

2m / 



that is easily solved numerically as an initial value problem [18]. 



However, in the most part of the NRQCD lattice evaluations only mean field contributions 



(devised to correct for the large tadpole terms in lattice perturbation theory [ 19 1) are considered 
while for the rest the matching coefficients are taken at the tree level. Such a procedure appears 
to be inconsistent: once the power counting in v is established order 0(a s ) corrections to the 
matching coefficients have to be included via at least one loop calculations. Notice also that, 



as discussed above, the loop corrections to the matching coefficients may be large [17] 



In NRQCD still the dominant role of the potential as well as the quantum mechanical nature of the 
problem (Schrodinger equation) are not maximally exploited: the NRQCD quark Green function 
satisfies in fact an equation (|l) in which potential and retardation contributions appear still mixed 



together |g, |2 



This goal is achieved by another EFT called potential Non-Relativistic QCD (pNRQCD) 



23, H 31] 



5 pNRQCD 

We want to build an effective field theory that describes the low energy region of the non-relativistic 
bound state, i.e. we want an EFT where only the ultrasoft degrees of freedom remain dynamical 
while all the other nonrelevant scales have been integrated out. This fixes two ultraviolet cutoffs 
for the new EFT: Ai and A2. The first one satisfies mv 2 <C Ai <C mv and is the cutoff of the 
energy of the quarks and of the energy and the momentum of the gluons, the second one satisfies 
mv <C A2 <C mv 2 and is the cutoff of the relative three-momentum p of the heavy quark bound 
system. This is the maximal reduction (=maximal number of degrees of freedom that we can make 
nondynamical) we can achieve in the description of the heavy quark bound state: in fact the heavy 
quarks with US energy have a soft three-momentum due to their non-relativistic dispersion relation. 

First we want to consider a situation where the matching to the new EFT is still perturbative, i.e. 
a situation in which the soft scale is still a perturbative scale: mv 3> Aqcd- The nonperturbative 
matching situation: mv ~ Aqcd> where the matching cannot be carried out via an expansion in a s 
is discussed in Sec. 12. Roughly speaking, we can say that the lowest excitations of quarkonium 
belong to the first situation while the excited states belong to the second situation. Indeed, the 
typical radius of the bound state is proportional to the inverse of the soft scale r ~ 1/mv and thus 
for the lowest states the condition mv <C Aqcd may be fulfilled [|J]. 



>C x w/q) 



NRQCD pNRQCD 

Figure 2: Matching procedure: off shell Green functions expanded in a s , 1/m and in the multipole 
expansion are equated in NRQCD and in pNRQCD. The difference is encoded in the potential 
matching coefficients that depend nonanalytically on the scale that has been integrated out, in this 
case r. 



6 pNRQCD for mv > A QC d 

We denote by R = (xi + X2)/2 the center of mass of the QQ system and by r = xi — X2 the relative 
distance. At the scale of the matching // (mv » // 3> mv 2 , Aqcd) we have still quarks and gluons. 
The effective degrees of freedom are: QQ states (that can be decomposed into a singlet <S(R, r, t) 
and an octet 0(R, r, t) under color transformations) with energy of order of the next relevant scale, 
0(Aqcd, rnv 2 ) and momentum p of order 0(mv), plus ultrasoft gluons Au(R, t) with energy and 
momentum of order O(AqcD) mv 2 ). Notice that all the gluon fields are multipole expanded. The 
Lagrangian is then an expansion in the small quantities p/m, 1/rm and in 0(Aqcd> mv 2 ) x r. 



6.1 The pNRQCD Lagrangian 

At the next-to-leading order (NLO) in the multipole expansion, i.e. at 0(r), we get |23|, 31 



L 



(i) 

pNRQCD 



1 F a F ^a + T] . jgf 



El 

m 



n=l 



m' 



(7) 



+Tr{ t (iDo-^-V -J2^-)o 



+gV A Tr {0 f r • E S + S f r • E o} + 5-yTr {o+r • E O + O f Or • e} . 

All the gauge fields in Eq. (^) are evaluated in R and t. In particular E = E(R, t) and iD$0 = 
id^O — <7[Aq(R, t),0]. The quantities denoted by Vj are the matching coefficients. They are functions 
of fi, fj/, r, p, Si, S2- We call V s and V Q the singlet and octet static matching potentials respectively. 
The equivalence of pNRQCD to NRQCD, and hence to QCD, is enforced by requiring the Green 
functions of both effective theories to be equal (matching). In practice, appropriate off shell ampli- 
tudes are compared in NRQCD and in pNRQCD, order by order in the expansion in 1/m, a s and 
in the multipole expansion. The difference is encoded in potential-like matching coefficients that 
depend non-analytically on the scale that has been integrated out (in this case r). 

At the leading order (LO) in the multipole expansion, the equations of motion of the singlet 
field is the Schrddinger equation 



i0 O 



+ V S 



s. 



m 



(8) 



Therefore pNRQCD has made explicit the dominant role of the potential and the quantum mechan- 
ical nature of the bound state. In particular both the kinetic energy and the potential count as 



mv 



in the v power counting. The leading orderf^] problem in fact reduces to the usual Schrddinger 
equation; the actual bound state calculation turns out to be very similar to a standard quantum 
mechanical calculation, the only difference being that the wave function field couples to US gluons 



3 The non-relativistic limit described by the Schrodinger equation with the static potential is called 'leading order' 
(LO); contributions corresponding to corrections of order v n to this limit are called N"LO. LL means 'leading log'. 



in a field theoretical fashion. From the solution of the Schrodinger equation come the leading order 
propagators for the singlet and the octet state, while the vertices come from the interaction terms 
at the NLO in the multipole expansion. In Fig.(||) you find such the Feynman rules in the static 
case (useful for the matching). For calculations inside pNRQCD the kinetic term p 2 /m has to be 
included in the singlet and octet propagators. 

The last line of (]?]) contains retardation (or nonpotential) effects that start at the NLO in the 
multipole expansion. At this order the nonpotential effects come from the singlet-octet interaction 
and the octet-octet interaction mediated by a ultrasoft chromoelectric field. 

Recalling that r ~ l/mv and that the operators count like the next relevant scale, 0(mv 2 , Aqcd)j 
to the power of the dimension, it follows that each term in the pNRQCD Lagrangian has a definite 
power counting. This feature makes L p nrqcd the most suitable tool for a bound state calculation: 
being interested in knowing the energy levels up to some power v n , we just need to evaluate the 
contributions of this size in the Lagrangian. 

We stress that 

• pNRQCD is equivalent to QCD. 

• pNRQCD has explicit potential terms and thus it embraces a description of heavy quarkonium 
in terms of potentials. 

• pNRQCD has explicit dynamical ultrasoft gluons and thus it describes nonpotential (retar- 
dation) effects. US gluons are incorporated in a second-quantized, gauge-invariant and sys- 
tematic way. In the appropriate dynamical situation Aqcd *C mv2 and Aqcd ~ mv 2 these 
give nonpotential terms of nonperturbative nature, like the term due to the gluon condensate 
||. From the power counting it follows that the interaction of quarks with ultrasoft gluons is 
suppressed in the lagrangian by v with respect to the LO (by gv if mv 2 ^> Aqcd)- 

• The power counting is unambiguous. Being all the scales disentangled, the perturbative 
calculation of the bound state is considerably simplified. 

• Calculations can be performed systematically in the v expansion and can be improved at the 
desired order. 

• Perturbative (high energy) and nonperturbative (low energy) contributions are disentangled, 
we can therefore systematically parameterize the nonperturbative contributions that we are 
not able to evaluate directly. 

• Lorentz invariance is again implemented through reparameterization invariance. In this case 
this implies the existence of relations between the potential matching coefficients [27]. 



• We have a definite prescription (power counting) on the calculations of the contribution of 
the various relativistic corrections to the energy levels. We have to solve first the Schrodinger 
equation with the static potential. All the rest, if it is suppressed by the power counting, has to 
be computed using quantum mechanical perturbation theory on the leading order Schrodinger 
solution. E.g. the spin-spin term proportional to the delta function I considered in Sec. 2 has 
to be calculated in quantum mechanical perturbation theory and has not to be included in the 
Schrodinger equation. 

In pNRQCD the potentials are defined upon integration of all the scales up to the ultrasoft scale 
mv 2 and are understood in the modern acception of matching coefficients: they are dependent on 
the scale of the matching // . Of course, such scale dependence is canceled in the energy levels by the 
contribution of the ultrasoft gluons which are cutoff at the same scale [)! . In particular pNRQCD 
provides us with a well defined procedure to obtaining the matching potentials via the matching at 
any order of the perturbative expansion. This is particularly relevant in QCD, where the calculation 



of the static potential is unclear in a naive perturbative frame, see e.g. [30, [H]]. Up to now the EFT 
approach is the only one able to supply a perturbative definition of the static potential at all orders, 
in particular to enable a concrete calculation at three loops (LL) and consequently a perturbative 
calculation of the energy levels at NNNLL. This is an example where the EFT seems not only useful 
but necessary. 
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Figure 3: Propagators and vertices of the pNRQCD Lagrangian (at order l/m° and at the NLO in 
the multipole expansion). 



6.2 The QQ potential in perturbation theory 

The existence of the different physical scales makes even a purely perturbative definition of the 
static QQ potential not free from complications. Let us consider the energy of the static quark 

sources E s (r) = lim — ln(Wn), (being Wo the static Wilson loop of size r x T, and the symbol 

T — >oo x 

( ) being the average over the gauge fields), which is usually considered as a definition of the static 
potential. At three loops E s shows infrared divergences |30|, |3l]| . These singularities may indeed be 
regulated, upon resummation of a certain class of diagrams, which give rise to a sort of dynamical 
cut-off provided by the difference between the singlet and the octet potential. However, such a 
dynamical scale is of the same order of the kinetic energy for quarks of large but finite mass and, 
therefore, should not be included into a proper definition of the static potential in the sense of the 
Schrodinger equation. This is similar to what happens for the Lamb shift in QED at order 1/m 2 . 
In QCD this effect calls, even in the definition of the static potential, for a rigorous treatment of 
the bound-state scales. To address the multiscale dynamics of the heavy quark bound state, the 
concept of effective field theory turns out to be not only helpful but actually necessary. 

We answer these questions by performing explicitly the singlet matching at order 1/mP and at 
the NLO in the multipole expansion. 

The matching can be done once the interpolating fields for S and O a have been identified 
in NRQCD. The former need to have the same quantum numbers and the same transformation 
properties as the latter. The correspondence is not one-to-one. Given an interpolating field in 
NRQCD, there is an infinite number of combinations of singlet and octet wave- functions with 
ultrasoft fields, which have the same quantum numbers and, therefore, have a non vanishing overlap 
with the NRQCD operator. However, the operators in pNRQCD can be organized according to the 
counting of the multipole expansion. For instance, for the singlet we have 



X f (x 2 , t)0(x 2 , xi ; t)V(xi, t) = Zy 2 (r)5(R, r, t) + Z l £ s (r)r r • E a (R, i)O a (R, r, t) + . . . , (9) 



and for the octet 

X f (x 2 , t)0(x 2 , R; *)r°0(R, xi ; i)V(xi, t) 



= Zy 2 (r)O a (K,r,t) 

+ Z^(r)rr-E a (R,t)5(R,r,t) + ..., 



(10) 



</>(y,x;t) = Pexp{ig J ds (y— x)-A(x— s(x— y), t)}, Zj being normalization factors. These operators 
guarantee a leading overlap with the singlet and the octet wave-functions respectively. Higher order 
corrections are suppressed in the multipole expansion. The expressions for the octet can be made 
manifestly gauge-invariant by inserting a chromomagnetic field in place of the pure color matrix. 
The fact that the matching can be done in a completely gauge-invariant way enables us to generalize 
pNRQCD to the case in which Aqcd — Tnv, i.e. to the nonperturbative matching, cf. Sec. 12. 

Now, in order to get the singlet potential, we compare 4-quark Green functions. From (^), we 
take in NRQCD 

/ = <)- 3 (x 1 -y 1 )5 3 (x 2 -y 2 )(^ n ), (11) 

and we equate (11) to the singlet propagator in pNRQCD at NLO in the multipole expansion (cf. 
Fig. 4 for a diagrammatic representation) 
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where <^> ad j is a Schwinger (straight-line) string in the adjoint representation and fields with only 
temporal argument are evaluated in the centre-of-mass coordinate. 
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NRQCD pNRQCD 

Figure 4: The matching between the Wilson loop in NRQCD and the singlet propagator in pNRQCD. 



Comparing Eqs. flllf ) and fll2f), one gets at the next-to-leading order in the multipole expansion 
the singlet wave- function normalization Z s and the singlet static potential V s . Va and V must have 
been previously obtained from the matching of suitable operators, but for the present purposes we 
only need the tree-level values: Va = 1 and V a = (Ca/2 — Cp)a s /r. Let us concentrate here on 
the matching potential V^. By substituting the chromoelectric field correlator in Eq. ( |i~2| ) with its 
perturbative expression we obtain at the next-to-leading order in the multipole expansion and at 
order at In a s 



V s (r) = £ s (r)| 2 _ loop+NNNLL + C F — — — In ^— . 



(13) 



The 2-loop contribution to E s has been calculated in 



. The NNNLL contributions arise from 
the diagrams studied first in |3C] and shown below. Note that V s and E s would coincide in QED 



C F C A< ln CAO 1 + 
12irr 2 



and that, therefore, this difference here is a genuine QCD feature. Such difference is switched on at 
NLO in the multipole expansion. An explicit calculation gives pl| 



V s (r) = -C F 



a v (r, //) 



(14) 
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where /3 n are the coefficients of the beta function (a s is in the MS scheme), and a\ and ct2 were 



given in [25]. We see that the interpretation of the potentials as matching coefficients in pNRQCD 
implies that the Coulomb potential is not simply coincident with the static energy E s . The Coulomb 
potential turns out to be sensitive to the ultrasoft scale. The same happens with the other potentials 
(like V or the potentials that bear corrections in l/m n ) that can equally be calculated via the 
matching procedure. 



6.3 The singlet potential in the situation mv 3> Aq CD ^> mv 2 

Since in this situation there is a physical scale (Aqcd) above the US scale, a potential can be properly 
defined only once this scale has been integrated out. At the NLO in the multipole expansion we get 



V s (r) = -C F 



ay(r, n') 



N r . 



M^T F Vl(r)- / dte-^-^)(E a (t)^(t,0)^E 6 (0))(M') 



(15) 



Therefore, the heavy quarkonium static potential V s is given in this situation by the sum of the 



purely perturbative piece calculated in Eq. (14) and a term carrying nonperturbative contributions 
(contained into non-local gluon field correlators). This last one can be organized as a series of power 
of r n by expanding exp{— it(V — V s )} = 1 — it(V — V s ) + . . . (since t ~ 1/Aqcd, V a — V s ~ mv 2 ). 
Typically the nonperturbative piece of Eq. dl5| ) absorbs the // dependence of ay so that the re- 
sulting potential V s is now scale independent. 

We notice that the leading nonperturbative term could be as important as the perturbative po- 
tential once the power counting is established and, if so, it should be kept exact when solving the 
Schrodinger equation. In Table 1 we summarize the different kinematic situations. 



mv 


mv' 2 


potential 


ultrasoft corrections 


» Aqcd 
» Aqcd 
» Aqcd 


> Aqcd 
~ Aqcd 
< Aqcd 


perturbative 
perturbative 
perturbative + 
short-range nonpert. 


US gluons + local condensates 
US gluons + non-local condensates 
No US (if light quarks 
are not considered) 



Table 1: Summary of the different kinematic situations. 



6.4 The perturbative potential and nonpotential nonperturbative corrections: 
Coulombic and quasi-Coulombic systems 

Let me define what I mean with Coulombic or quasi-Coulombic systems |37|. If mv S> mv 2 > Aqcdj 
the system is described up to order af by a potential which is entirely accessible to perturbative 
QCD (|l4|). Nonpotential effects start at order ln^' J32"[ . We call Coulombic this kind of system. 
Nonperturbative effects are of nonpotential type and can be encoded into local (a la Voloshin- 
Leutwyler(B|) condensates (if mv 2 S> Aqcd) or non-local condensates (if mv 2 ~ Aqcd): they are 
suppressed by powers of Aqcd /mv 2 and Aqcd/"^ respectively. We will discuss in details the case 
of Coulombic systems in the next section. 

If mv S> Aqcd ^ mv 2 , the scale mv can be still integrated out perturbatively, giving rise 
to the Coulomb- type potential fli~4|), Nonperturbative contributions to the potential arise when 



integrating out the scale Aqcd [23], precisely as explained in Sec. 6.3. We call quasi- Coulombic the 
systems where the nonperturbative piece of the potential can be considered small with respect to 
the Coulombic one and treated as a perturbation. 

Some levels of tt, the lowest level of bb may be considered Coulombic systems while the 
J/tp, the n c and the short-range hybrids may be considered quasi-Coulombic. The B c may be in a 
boundary situation. As it is typical in an effective theory, only the actual calculation may confirm if 
the initial assumption about the physical system was appropriate. In Sec. 9 I will give an example 
of computation in a quasi-Coulombic system after an appropriate redefinition of the mass parameter 
that subtracts the largest indetermination in the perturbative series. 

For all these systems it is relevant to obtain a determination of the energy levels as accurate as 
possible in perturbation theory. 



7 CALCULATION of QUARKONIUM ENERGIES at 0{ma 5 s In a s ) 

The perturbative energy levels of quarkonium are known at 0(maf)[^]. pNRQCD provides us with 
a well defined way of calculating the energy levels at higher orders, since the size of each term in 
the pNRQCD Lagrangian is well-defined and since we are now able to calculate the static potential 
beyond two loops. In order to obtain the leading logs at 0{ma^) in the spectrum, V s has to be 
computed at 0(a^ In), vj 1 ^ at 0(0^ In), uj 2 ** at 0(a 2 In) and vj 3 ^ at O(o s ln). The matching to 
pNRQCD at N 3 LL accuracy was calculated in [^] and the potentials at the requested accuracy 
were thus obtained. 

The total correction to the energy at 0(maf In a s ) is given by the sum of the averaged values 
of the potentials plus the nonpotential contributions, 

SE nM = 5 pot E n ij(fi>) + * us £„,,(m') , (16) 




In 



A* 



ma* 



(17) 



where E n = -mC 2 F a 2 J {An 2 ) and G jjL = -(I + 1)(4/ - 1)/(2Z - 1) for j = I - 1, = -lfor j = I and 
= /(4Z + 5)/(2Z + 3)forj = 1 + 1. 

The lna s appearing in Eq. (|l7| ) come from logs of the type lnl/mr in the potential. Therefore 
they can be deduced once the dependence on lnm is known. The \j! dependence of Eq. (|T^) cancels 
against contributions from US energies. At the next-to-leading order in the multipole expansion 
the contribution from these scales reads 



S VS E n ^') 



3Nr. 



poo 

Tpx / dt{nJ\re it{E "- H ° ) r\nJ)(E a (t)<P(t,0)^ 3 E\0))(u') 



(18) 



where H Q = p 2 /2m + V a . 

Different possibilities appear depending on the relative size of Aqcd with respect to the US 



scale maj . If we consider that Aqcd 



~ ma' 



the gluonic correlator in Eq. (^) cannot be computed 



using perturbation theory. We are still able to obtain all the ma\ ln(ma s /m) and ma\ \n(ma s / u') 
contributions where the \J dependence cancels now against the US contributions that have to be 
evaluated non-perturbatively. 



Actually both bb and cc ground states have been studied in this way|33 



can be computed perturbatively. Being ma 2 the 



If we consider that ma s 2 S> Aqcd> Eq. (f 
next relevant scale, the effective role of Eq. (|18| ) will be to replace fj! by ma 2 (up to finite pieces 
that we are neglecting) in Eq. Then Eq. ( |l6| ) simplifies to 
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(19) 



plus non-perturbative corrections that can be parameterized by local condensates |32| , |5[ and are 
of order a 2 (^^) 2 (^£) 2 and higher. For the T(15) and tt non-perturbative contributions are 
expected not to exceed 100 150 MeV and 10 MeV respectively. 

The calculation (|l^) paves the way to the full N 3 LO order analysis of Coulombic systems and 
is relevant at least for tt production and T physics. In the first case it is a step forward reaching a 
50 MeV sensitivity on the top quark mass for the ti cross section near threshold to be measured at 
future Linear Colliders |]43[| . In the second case they improve our knowledge of the b mass |33| . 

In particular, from ( |l9| ) we can estimate the N 3 LL correction to the energy level of the T(1S) 
and we find [32] 
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m b a 4 s (n)a s (n') ]nl/a a (pf) ~ (80 4- 100) MeV, 



(20) 



which appears not to be small. Corrections from this N 3 LL terms have been calculated for the tt 



cross section and T(1S) wave function, cf. |35j and |33j. They also turn out to be sizeable. 

Large corrections, however, already show up at NLO and N 2 LO and are responsible for the bad 
convergence of the perturbative series in terms of the pole mass. 



8 RENORMALONS, the POLE MASS and the PERTURBATIVE 
EXPANSION 



The bad convergence of the perturbative expansion can be, at least in part, attributed to renormalon 
contributions. The pole mass, thought an infrared safe quantity^], has long distance contributions 
of order Aqqd [flCf] . Also the static potential is affected by renormalons (see e.g. [@). Rephrasing 
them in the effective field theory language of pNRQCD we can say that the Coulomb potential 
suffers from IR renormalons ambiguities with the following structure 



V s (r)\ 1Rmn = C + C 2 r 2 + . . . 



(21) 
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Figure 5: Renormalon in the static potential at LO in the multipole expansion, here \i is a hard 
cut-off, 1/r > n > Aqcd- 
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Figure 6: Renormalon contribution the pole mass. 

The constant Cq ~ Aqcd is known to be canceled by the IR pole mass renormalon (2m po i e |iR ren = 
—Co, (4^] ). Several mass definitions appropriate to explicitly realize this renormalon cancellation 
have been proposed^]. Among the others, the IS mass|3^] is defined as half of the perturbative 
contribution to the 3 Si qq mass. Unlike the pole mass, the IS mass, containing, by construction, 
half of the total static energy (2m + y Coul ^ [ s free Q f ambiguities of order Aqcd- Taking e.g. the 
T(1S), the IS" mass is related to the physical T(1S) mass by E(T(1S)) = 2mis + Ax- Ax is the 
poorly known non-perturbative contribution, which is likely, as we said, to be less than 100 -j- 150 
MeV. In the next section we will present an explicit example that shows how, using this mass and 
thus dealing with quantities that are infrared safe at order Aqcd , the pathologies of the perturbative 
series, due to the renormalon ambiguities affecting the pole mass, are cured. 

It is possible to show that the second infrared renormalon, C2 — Aq CD , of V s cancels against 
the appropriate pNRQCD UV renormalon in the contribution to the potential originating at NLO 
in the multipole expansion. What remains is the explicit expression for the operator which absorbs 
the C2 ~ Aq CD ambiguity (for details see f23|): such nonperturbative operator turns out to be a 
nonlocal chromoelectric condensate of the type used in some QCD vacuum models [37, [H| ]. 



An interesting open question is if an explicit renormalon subtraction similar to that one at 
O(Aqcd) can be realized at the subsequent order 0(Aq CD ). This may be relevant since this renor- 
malon is related to the corrections at N 3 LL discussed in the previous section. Indeed, it is still 
an open problem whether the largeness of the N 3 LL corrections is an artifact due to our partial 
knowledge of the contributions at this order, or if it is an artifact due to the fact that the sub- 
sequent renormalon cancellation has to be realized at this order or finally if it is a true signal of 
the breakdown of the perturbative series. To make more definite statements one should know the 



complete N 3 LO or understand the mechanism of cancellation of the second renormalon [42|. 

In the next section I will present a concrete example of the relevance of the mass renormalon 
cancellation in order to obtain reliable phenomenological predictions. 

9 A QUASI COULOMBIC SYSTEM with 

EXPLICIT O(Aqcd) RENORMALON CANCELLATION 

I consider here the perturbative calculation up to order ma\ of the energy of the be ground state: 
this will be relevant to a QCD determination of the B c mass if this system is Coulombic or at least 
quasi-Coulombic. I will assume this to be the case. 

In order to calculate the B c mass in perturbation theory up to order af, we need to consider 
the following contributions to the potential: the perturbative static potential at two loops, the 1/m 
relativistic corrections at one loop, the spin-independent 1/m 2 relativistic corrections at tree level 
and the 1 / m 3 correction to the kinetic energy. We do not consider In a s corrections since the 



mechanism responsible for this large contributions has not yet been understood. Then, we have 1 46] 



E(B c ) pert =m b + m c + E (fi) jl- ^ 
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(22) 



being m re d = mhm c / '(mj + m c ) , Z = ln(2Ci?a s m re d//i)) Eo(p) = — m re d(C'i7'a s (^)) 2 /2 and /2 the scale 
around which we expand a s (r). 

If we use here the pole masses m& = 5 GeV, m c = 1.8 GeV and p = 1.6 GeV, then we obtain 



E(B c ) pert ~ 6149 MeV ~ 6800 - 115 - 183 - 353 MeV, where the second, third and fourth figures 
are the corrections of order a^, and af respectively. The series turns out to be very badly 
convergent. This reflects also in a strong dependence on the normalization scale p: at p = 1.2 GeV 
we would get E(B c ) pert ~ 5860 MeV, while at p = 2.0 GeV we would get E(B c ) pcrt ~ 6279 MeV. 
The non-convergence of the perturbative series (|22|) signals the fact that large /3q contributions 
(coming from the static potential renormalon) are not summed up and canceled against the pole 
masses. In order to obtain a well-behaved perturbative expansion, we use, now, the so-called IS" 
mass. We consider the perturbative contribution (up to order ctf ) of the 3 Si levels of charmonium 
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Figure 7: E(B c ) pcj - t as a function of fi at A^ 
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and bottomonium: 



E{J/^) P ert = f(m c ); £(T(lS)) pert = f(m b ) 



which are respectively a function of the c and the b pole mass and can be read off from Eq. ( |22| ) 
in the equal-mass case, adding to it the spin-spin interaction energy: m{Cpa^) /3. We invert these 
relations in order to obtain the pole masses as a formal perturbative expansion depending on the 
IS masses. Finally, we insert the expressions m c = f^ 1 (E(J/ip) per t) and m& = f~ 1 (E(T(lS)) per t) 
in Eq. (^). At this point we have the perturbative mass of the B c as a function of the J /if) and 
T(15) perturbative masses. If we identify the perturbative masses E(J/ip) pev t, E(T(lS)) peT t with 
the physical ones, then the expansion (^) depends only on the scale p. The perturbative series 
turns out to be reliable for values of p bigger than (1.2 -j- 1.3) GeV and lower than (2.6 -f- 2.8) 
GeV. For instance, E(B c ) pevt = 6278.5 + 35 + 6.5 + 5.5 MeV at the scale p = 1.6 GeV. Therefore, 
we obtain now a better convergence of the perturbative expansion and a stable determination of 
the perturbative mass of the B c . This fact seems to support the B c being indeed a Coulombic or 
quasi-Coulombic system. By varying p from 1.2 GeV to 2.0 GeV and A^ from 250 MeV to 350 



MeV and by calculating the maximum variation of E(B C 
get as our final result 



E(B C ) 



port 



6326±g 9 MeV 



ip er t in the given range of parameters, we 

(23) 



As a consequence of the now obtained good behaviour of the perturbative series in the considered 
range of parameters, the result appears stable with respect to variations of p (see Fig. 0) and, 
therefore, reliable from the perturbative point of view. It represents a rather clean prediction of the 
lowest mass of the B c . Notice that all the existing predictions are based either on potential models 
or on lattice evaluation (with still large errors). 

Non-perturbative contributions have not been taken into account so far. They affect the iden- 
tification of the perturbative masses E(B c ) pert , E(T(lS)) peT t, E(J/ip) pert , with the corresponding 
physical ones through Eq.(22). Let us call these non-perturbative contributions A# c , Ax and Aj/^ 



respectively. As discussed before, they can be of potential or nonpotential nature. In the latter 
case they can be encoded into non-local condensates or into local condensates. Non-perturbative 
contributions affect the identification with the physical B c mass roughly by an amount ~ —Aj/,^/2 
-At/2 +A Bc . Assuming < 300 MeV, |A T | < 100 MeV and A T < A Bc < A J/V ,, the identifi- 

cation of our result ( f23|) with the physical B c mass may, in principle, be affected by uncertainties, 
due to the unknown non-perturbative contributions, as big as ±200 MeV. However, the different 
non-perturbative contributions A are correlated, so that we expect, indeed, smaller uncertainties. 
If we assume, for instance, At and A.j/^ to have the same sign, which seems to be quite reasonable, 
then the above uncertainty reduces to ±100 MeV. This would confirm, indeed, that the effect of 
the non-perturbative contributions on the result of Eq. ( |2"3| ) is not too large. 



10 RENORMALIZATION GROUP IMPROVEMENT 
and vNRQCD 

Of course the results presented here can be renormalization group improved, i.e. the potentially 
large logarithms can be resummed. The RG group improvement in pNRQCD has been done only at 
the level of the static potential E^] up to now. When going beyond the static limit v = a special 
care has to be payed to the correct implementation of the renormalization group equations because 
on the bound state we are dealing with correlated scales mv and mv 2 Another EFT, called 



vNRQCD [26], has been devised explicitly to deal with the resummation of large logs. In this case 
the treatment is, however, purely perturbative. The Lagrangian is explicitly decomposed in soft, 
potential and ultrasoft fields: 

L V NRQCD = L s + Lp + L u 

The matching is performed at the scale m. The RG equations evolve simultaneously from the scale 
m to the correlated scales mv and mv 2 . This seems to provide the correct double logarithms. 

The approach of vNRQCD may be complicated by the fact that the different scales remain 
entangled. Moreover vNRQCD is perturbative by construction and thus cannot be used in order 
to parameterize nonperturbative effects. The goal of vNRQCD is summing up the logarithms 
in the velocity. In my opinion, however, the real challenge in heavy quark systems remain the 
nonperturbative corrections. It is also important to keep in mind that, as I noticed in Sec. 4, the 
finite part contributions may be as much (or more) relevant as the logarithmic ones. 



11 RESULTS and APPLICATIONS 

In this Section we just recall a list of physical situations where the presented results of pNRQCD 
(in the situation of the perturbative matching mv 3> Aqcd) have been or may be relevant: 

• Toponium cross section J^3| , 

• Determination of the mass of the b, once the renormalon contribution has been appropriately 
subtracted g|, E§> 0- 

• Hybrids and gluelumps 

pNRQCD gives model independent predictions on the behaviour of the hybrid static potentials. 
In particular it predicts for these potentials an octet behaviour at very short distances and it 
correctly states all the degeneration patterns in the small r region (2^, |45j. Moreover, it allows 
to relate the mass of the gluelumps to the correlation lengths of some nonlocal vacuum field 
strength correlators |37], |23| allowing us to obtain interesting model independent information 
on the behaviour of these nonperturbative objects |^4j . 



Quarkonium scattering, Van Der Waals forces, Quarkonium Production. Work is in progress 
on these applications. 



12 NONPERTURBATIVE MATCHING 
pNRQCD for A~mw 

In this case the potential interaction is dominated by nonperturbative effects. This is the most 
interesting situation, in which most of the mesons seem to lie. 

A large effort has been made in the last decades in order to obtain from QCD the nonpertur- 
bative potentials in the Wilson loop approach. pNRQCD allows us to obtain via a nonperturbative 
matching all the nonperturbative potentials |j47| . On this respect I want to discuss few concepts 
and results. 

In pNRQCD a potential picture for heavy quarkonium emerges at the leading order in the US 
expansion under the condition that the matching between NRQCD and pNRQCD can be performed 
within an expansion in 1/m. The gluonic excitations (hybrids and glueballs) that form a gap of 
order Aqcd with respect to the quarkonium state can be integrated out and the potentials follow 
in an unambiguous and systematic way from the nonperturbative matching to pNRQCD. Thus, we 



recover the quark model from pNRQCD [47]. The US degrees of freedom in this case are not coloured 
gluons but US gluonic excitations between heavy quarks and pions. They can be systematically 
included and may eventually affect the leading potential picture. Let us consider for instance the 
singlet matching potential. Disregarding the US corrections, we have the identification 

V s (r) = lim ^\n{Wo). 

US corrections to this formula are due to pions and US gluonic excitations between heavy quarks, 
and may be included in the same way as the effects due to US gluons have been included in the 
perturbative situation. 

The complete 1/m 2 potential has been calculated along these lines [47]. There are many appeal- 
ing and interesting features of this procedure. The matching calculation is performed in the way 
of the quantum mechanical perturbation theory on the QCD Hamiltonian (where perturbations are 
counted in orders of 1/m and not of a s ) and only at a later stage the relation with the Wilson 
loop and field insertions is established. This allows us to have a control on the Fock states of the 
problem and on the contributions coming from gluonic excitations in the intermediate states. At 
the end, all the expressions are again given in terms of Wilson loops, which can be evaluated on 
the lattice or in QCD vacuum models. The potentials turn out to be naturally factorized in a hard 
part (the matching coefficients at the hard scales inherited by pNRQCD from NRQCD) and a low 
energy part (the Wilson loops expressions). The power counting may turn out to be quite different 
from the perturbative (QED-like) situation and this may turn out to be very important for some 
applications (e.g. quarkonium production). 

I will not discuss these very interesting results any longer because it is the content of the review 
of Antonio Vairo at this Conference [47]. 



13 CONCLUSION and OUTLOOK 

I have shown that it is possible to construct systematically and within a controlled expansion an 
effective theory of QCD, which describes heavy quark bound states. All known perturbative and 
nonperturbative regimes (potential, nonpotential) , are dynamically present in the theory, which is 
equivalent to QCD. I have presented many applications of pNRQCD in the situation Aqcd *C mv - 
In the situation Aqcd — mv, I have shortly discussed how pNRQCD allows us to systematically 



factorize the nonperturbative heavy quark dynamics. Similar techniques have been applied to 
pionium |4S[ |. 
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